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ABSTRACT 

Context. It is still debated whether or not the Galactic chemical evolution of sulphur in the halo followed the constant or flat trend 
with [Fe/H], ascribed to the result of explosive nucleosynthesis in type II SNe. It has been suggested that the disagreement between 
different investigations of sulphur abundances in halo stars might be due to problems with the diagnostics used, that a new production 
source of sulphur might be needed in the early Universe, like hypernovae, or that the deposition of supernova ejecta into the interstellar 
medium is time-delayed. 

Aims. The aim of this study is to try to clarify this situation by measuring the sulphur abundance in a sample of halo giants using 
two diagnostics; the S i triplet around 1045 nm and the [S i] line at 1082 nm. The latter of the two is not believed to be sensitive to 
non-LTE effects. We can thereby minimize the uncertainties in the diagnostic used and estimate the usefulness of the triplet in sulphur 
determination in halo K giants. We will also be able to compare our sulphur abundance differences from the two diagnostics with the 
expected non-LTE effects in the 1045 nm triplet previously calculated by others. 

Methods. High-resolution near-infrared spectra of ten K giants were recorded using the spectrometer CRIRES mounted on VLT. Two 
standard settings were used; one covering the S i triplet and one covering the [S i] line. The sulphur abundances were determined 
individually with equivalent widths and synthetic spectra for the two diagnostics using tailored ID model atmospheres and relying on 
non-LTE corrections from the litterature. Effects of convective inhomogeneities in the stellar atmospheres are investigated. 
Results. The sulphur abundances derived from both the 1082 nm [S i] line and the non-LTE corrected 1045 nm triplet favor a flat 
trend for the evolution of sulphur. No 'high' values of [S/Fe] shown in some previous studies are seen in our sample. 
Conclusions. We corroborate the flat trend in the [S/Fe] vs. [Fe/H] plot for halo stars found in other works and cannot find a scatter 
nor a rise in [S/Fe] obtained in some other previous studies. We find the sulphur abundances deduced from the non-LTE corrected 
triplet somewhat lower than the abundances from the [S i] line, possibly indicating too large non-LTE corrections. Considering 3D 
modeling, however, they might instead be too small. Further we show that the [S i] line is possible to use as a sulphur diagnostic down 
to [Fe/H] 2.3 in giants. 
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1. Introduction 

Sulphur is an ff-element (like O, Ne, Mg, Si, Ar, and Ca) and 
these are believed to be produced mainly in SNe type II by addi- 
tions of a-particles. The SNe type la, on the other hand, mainly 
produce the iron peak elements (Fe, Co, and Ni). Due to the 
different life times of the two groups of productions sites, the in- 
terstellar medium (ISM) in the early Galaxy is believed to con- 
tain higher abundance a-elements compared to Fe than later on 
when the SNe type la start expelling iron. Thus, stars that form 
at a certain time serve as markers of the [a/Fe J] ratio in the 
ISM at that particular stage in the evolution of the Galaxy. This 
requires that no extra a-elements nor Fe has been produced in 
the star and contaminated the photosphere during the star's life. 
The theory of nucleosynthesis and stellar structure and evolution 
suggest this to be the case. Since the iron abundance, [Fe/H], 



* Based on observations collected at the European Southern 
Observatory, Chile (ESO program 080.D-0675(A)). 

1 The notation [A/B] = log (N A /N B ), - log (N A /N B ) Q , where N A and 
Nb are the number abundances of elements A and B respectively. 



roughly increases with time, an [a/Fe] vs. [Fe/H] plot is ex- 
pected to show a plateau for the lowest metallicities and a neg- 
ative slope from the time when the SNe type la start producing 
large amounts of iron. Regarding sulphur in particular, it is be- 
lieved to be produced via oxygen burning just like Si and Ca 
and therefore these three elements are expected to vary in lock- 
step with each other. The expected flat behavior in the [a/Fe] vs. 
[Fe/H] plot is however not unambiguously observed for sulphur 
in all previous works and hence there is no agreement regard- 
ing its evolution. This mismatch with the rest of the a-elements 
makes the behavior of sulphur interesting to investigate for its 
own sake, but the Galactic evolution of sulphur is also important 
to determine for at least two other reasons: 

- The Galactic evolution of sulphur can be matched against the 
predictions from modeling of supernovae yields and thereby 
testing such models. 

- The volatile properties of sulphur mean that the gas abun- 
dances of sulphur measured in, for example, damped Lyman 
a-systems, believed to be galaxies in making, do not need to 
be corrected for dust depletion to estimate the sulphur abun- 
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dance in the system. This abundance together with Zn might 
be used as a proxy for the [a/Fe] and thereby the timescale 
of the star formation history of the system can be inferred 
(Nissen et al. 2004 ). In order to make such conclusions possi- 
ble the Galactic evolution of sulphur needs to be determined 
first. 

Despite its importance, the Galactic evolution of sulphur is 
not a well studied subject (compared to, e.g., the Galactic evolu- 
tion of oxygen). This is mainly due to the lack of suitable sulphur 
diagnostics. During the years several diagnostics have been used 
with different strengths and weaknesses. Sulphur, being situated 
directly below oxygen in the periodic table, has a similar set of 
energy levels and transitions as oxygen. For example the 1082 
nm [S i] line and the 1045 nm triplet used in this work are anal- 
ogous to the widely used [O i] line at 630 nm and the O i 845 nm 
lines. 

It turns out that different works, using different sulphur di- 
agnostics, have resulted in several different scenarios for the 
Galactic evolution of sulphur. From measurements of the 869 
nm S i doublet Israelian & Rebolo (2001 ) and Takada-Hidai et al.] 
(2002 ) found a negative slope in the [S/Fe] vs. [Fe/H] plot for all 
observed iron abundances (-3 < [Fe/H] < +0.5). To explain this 
they proposed either a scenario for sulphur production involving 
hypernovae (e.g., Nakamura et al. 2001) or time-delayed iron 
deposition into the ISM as compared to sulphur (see Ramaty 
2000} for a description of the mechanism in connection to 



etal. 



[O/Fe] vs. [Fe/H]). A majority of the stars observed by Israelian 
|& Rebolo| ( |200TT ) and |Takada-Hidai et al.| ( |200"2"l > were later re- 
analyzed using the roughly ten times stronger S i triplet around 
923 nm ( |Ryde & Lambert|2004||Korn & Rydep 005), allowing 
a more precise abundance determination for metal-poor stars. 
These new results, contradicting the earlier findings, suggest an 
evolution analogous to the rest of the a-elements. Neither diag- 
nostic is, however, optimal; the 923 nm triplet is stronger than 
the 869 nm doublet, but on the other hand it is situated in a spec- 
tral region heavily plagued with telluric lines and the middle line 
of the triplet is situated in the wing of the strong Paschen ( line. 

Further |Takada-Hidai et aL| ( |2005[ ) analyzed 22 stars using 
both the 923 nm triplet and the 869 nm doublet finding a flat 
trend from both diagnostics, but significant differences in abun- 
dances for the stars with [Fe/H] < -1.5. These differences 
might be due to non-LTE effects in the diagnostics used, and 



this was investigated by Takeda et al. (2005) by calculating non- 
LTE corrections for the 869 nm doublet, the 923 nm triplet, and 
in addition the 1045 nm triplet. They used data for determin 
ing the evolution of sulphur from Israelian & Rebolo (|2001 



Takada-Hi dai et al.| (|2002l>; |Ryde & Lambert| ( |2004| ); |Nissen 
et al.| ( |2004| i; |Takada-Hidai et al.| ( |2005| l and some older works 



and applied their non-LTE corrections to the 869 nm doublet and 
the 923 nm triplet (the 1045 nm triplet not being observed at the 
time). Surprisingly, they found that the discrepancies between 
the evolution of sulphur measured using the different diagnos- 
tics increased; the 869 nm doublet suggesting a steady increase 
of [S/Fe] for lower [Fe/H] and the 923 nm triplet indicating a 
plateau for halo stars. To resolve this mismatch and check their 
non-LTE modeling they proposed observations of the 1045 nm 
triplet. 



Caffau et al. ( 2005 i combined spectra from four observation 



runs, initially aimed at studying other elements into the largest 
sample of Galactic sulphur abundance measurements with a to- 
tal of 74 dwarfs. They used the NTT and VLT telescopes and 
since the wavelength coverage of the observation runs varied, 
they did not use the same diagnostics in all their sulphur abun- 



dance determinations. For each star they used as many diagnos- 
tics as possible of a very weak multiplet around 675 nm, the 
doublet around 869 nm, and the triplet around 923 nm. They 
found an [S/Fe] vs. [Fe/H] plot with a large scatter in [S/Fe] for 
-2.4 < [Fe/H] < — 1 resembling a combination of the two types 
of evolution previously described. 

The development of better infrared spectrometers has made 
more sulphur lines measurable and the problem with finding a 
suitable diagnostic for determining sulphur abundance is nowa- 
days less severe. One example is the [S i] line at 1082 nm used 
in this work. This line is not believed to be affected by non-LTE 
effects, but unfortunately the [S i] line is undetectable in halo 
dwarfs for low [Fe/H]. It is however detectable in giants down 
to [Fe/H] ~ -2.5. Another example is the triplet around 1045 
nm also used in this work. This triplet is not as strong as the 
triplet at 923 nm, but it is situated in a spectral region almost 
unaffected by telluric lines. Caffa u et al.| ( |2007] > investigated the 
possibility of using the 1045 nm triplet for sulphur determination 
in disk dwarfs and it was recently used by them for determining 
sulphur abundances in four halo dwarfs corroborating the scatter 
in [S/Fe] fo r [Fe/H] ~ -1 fou nd in |Caffau et aIJ ( |2"0T)5] l but at a 
lo wer level (|Caffau et al.|201Q|>. 

Takeda & Tak ada-Hidai| ( |20lT l also recently used the 1045 



nm triplet to determine [S/Fe] in 33 halo/disk stars proposing yet 
a new scenario for the evolution of sulphur; a zig-zag trend with 
a local plateau around [S/Fe] ~ 0.3 for -2.5 < [Fe/H] < -1.5 
preceded by a rise in [S/Fe] for lower [Fe/H]. This is however 
not confirmed by Spite et al. ( 201 1\ using the 923 nm triplet to 
determine sulphur abundance in 33 giants and turnoff stars, with 
iron abundances [Fe/H] < -2.5, without finding any stars with 
'high' [S/Fe]. 

Obviously there is no consensus on the Galactic chemical 
evolution of sulphur. Is there a scatter or some kind of rise for 
halo stars or not? In this work we will present and compare sul- 
phur abundance measurements for 10 halo K giants using both 
the [S i] line at 1082 nm and the 1045 nm triplet. Since the [S i] 
line is not affected by non-LTE effects we will also be able to 
estimate the non-LTE effects of the 1045 nm triplet by compar- 
ing the derived sulphur abundances from the two diagnostics. 
Thereby we will be able to check the validity of the non-LTE 



corrections calculated by Takeda et al. ( 2005[ ) for halo K giants. 
Since non-LTE corrections seem to be significant for most sul- 
phur diagnostics this empirical test of the corrections is impor- 
tant. 



2. Observations 

We have observed ten K giants in the Galactic halo using the 
spectrometer C RIRES (Ka ufl et al] [20041 |Moorwood| |20"0"5j 
Kaufl et al. 2006) mounted on VLT. Giants were chosen because 
the 1082 nm [S i] line is too weak to be observed in metal-poor 
dwarfs (see Sect. |2.1| i. CRIRES is a high-resolution echelle spec- 
trometer designed for near-infrared observations and it uses nod- 
ding and jittering to eliminate the sky background and adaptive 
optics to enhance the S/N. Basic data for the observed stars are 
shown in Table [T] and a summary of the observations is shown 
in Table [2] The difference in S/N in the two settings is in some 
cases due to different integration times and in changes due to 
clouds in the star's visibility between the two observations. 

The observations were made with a slit width of 0.25" re- 
sulting in a spectral resolution of R = A/AA ~ 80000 and 2.5 
pixels per resolution element. They were carried out in service 
mode during October 2007 - February 2008. We used two stan- 
dard settings: one covering the Si triplet around 1045 nm and 
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Table 1. Basic data for the observed stars. 



Star identifier 


RA (J2000) 1 


Dec (J2000) 1 


i j 1 

V 


B-V 


Parallax 1 






(h m s) 


(d am as) 








(mas) 


HD 13979 


HIP 10497 


02 15 20.8536 


-25 54 54.861 


9.17 


0.647 ± 


0.024 


1.93 ± 


1.25 


HD21581 


HIP16214 


03 28 54.4853 


-00 25 03.117 


8.70 


0.786 ± 


0.018 


4.27 + 


1.20 


HD23798 


HIP17639 


03 46 45.7217 


-30 51 13.329 


8.28 


1.033 + 


0.015 


0.71 ± 


0.95 


HD26297 


HIP19378 


04 09 03.4175 


-15 53 27.068 


7.46 


1.088 ± 


0.011 


1.28 ± 


1.01 


HD29574 


HIP21648 


04 38 55.7328 


-13 20 48.138 


8.33 


1.304 ± 


0.019 


0.66 ± 


1.03 


HD36702 


HIP25916 


05 31 52.2305 


-38 33 24.046 


8.33 


1.148 ± 


0.018 


0.68 ± 


0.78 


HD44007 


HIP29992 


06 18 48.5269 


-14 50 43.424 


8.05 


0.829 ± 


0.012 


5.17 ± 


1.02 


HD83212 


HIP47139 


09 36 19.9533 


-20 53 14.759 


8.33 


1.015 ± 


0.015 


1.96 ± 


0.98 


HD85773 


HIP48516 


09 53 39.2415 


-22 50 08.425 


9.42 


1.078 ± 


0.046 


4.07 ± 


1.30 


HD103545 


HIP58139 


11 55 27.1618 


+09 07 45.028 


9.42 


0.844 ± 


0.043 


0.60 ± 


1.23 



References. (l) |Perryman et al.|fl997) 



Table 2. Summary of the observations with CRIRES/VLT. 



Star 


Date 


Integration time" 


S/N* 








Triplet 


[Si] 


Triplet 


[Si] 






(s) 


(s) 






HD 13979 


30 Oct. 2007 


1080 


1440 


360 


570 


HD21581 


16 Jan. 2008 


960 


960 


330 


310 


HD23798 


31 Oct. 2007 


840 


840 


500 


520 


HD26297 


16 Jan. 2008 


240 


240 


200 


330 


HD29574 


16 Jan. 2008 


360 


360 


300 


310 


HD36702 


16 Jan. 2008 


720 


720 


320 


360 


HD44007 


15 Feb. 2008 


480 


480 


370 


350 


HD83212 


16 Jan. 2008 


480 


480 


310 


240 


HD85773 


16 Jan. 2008 


1080 


1080 


480 


440 


HD103545 


16 Jan. 2008 


1440 


1080 


110 


160 



Notes. 

" The total integration times are given by NDIT x DIT x 

NEXP x NABCYC LES x 2, see the CRIRES U ser Manual at 

http://www.eso.org/sci/facilities/paranal/instruments/crires/doc/ 

" Approximate signal-to-noise ratios per pixel of the observed spectra. 



one covering the [S i] line at 1082 nm. The spectral range of the 
two settings used are roughly 12 nm with the sulphur lines as 
centered as possible. 

A couple of fast rotating B stars were also observed used for 
checking that no telluric lines were affecting the lines analyzed. 

2.1. Blending lines and stellar sample 

Using only a few lines, or even just one, for an abundance de- 
termination such as ours calls for an extra careful examination 
of blends, and in case blends are found: to try to choose the 
stellar sample in such a way that the effects of the blending 
are minimized. The [O i] line at 630 nm, analogous to our [S i] 
line, is for example blended with a Ni i line, but the blending 
is not relevant for metal-poor stars ( |Allende Prieto et al.|2Q01) >. 
To determine for which of our sulphur lines there might be sig- 
nificant blends we have calculated synthetic equivalent widths 
for all lines in the relevant wavelength region in a grid of 288 
model atmospheres. For metals we have used a line list of a rele- 



vant wavelength section from the VALD I database ( Kupka et al. 
|1999| l updated according to |Gustafsson et al.| ( |2008| l. The lists 
of molecules we were provided by Bengt Edvardsson (private 
communication), in turn mostly compiled by Bertrand Plez, and 
they include CH ( |j0rgensen et al.||1996| >, OH ( |Goldman et al. 



1998|>, CrH (|Burrows et al.|2002|l, SiH (electro nic K urucffi, FeH 
(jDulick et al.|2003| l, H 2 Q ( [Barber et al.|2006| >, C 2 ( |Querci et al.] 
|1971[|1974| i, CaH, CN, TiO, VO, and, ZrO (electronic and un- 
published Ple^. The list of lines with equivalent widths greater 
than zero within ±0.05 nm (Av ~ 30 km/s) from the used sul- 
phur lines is presented in Table [3] The equivalent widths of the 
1082 nm [Si], the 1045.5 nm and 1045.7 nm triplet lines, and 
blends with strengths in parity with the relevant sulphur line are 
shown in Figs. [T]- [3] The 1045.9 nm triplet line does not have 
any known signicant blends. 

The only blend for the [S i] line present in halo giants is that 
due to Cr i, but this line is separated enough from the sulphur line 
to be resolvable in our high resolution spectra, see for example 
the Cr i line in Fig. [5] For the 1045.5 nm and 1045.7 nm S i triplet 
lines the 1045.54 nm and 1045.69 nm Fe i lines shown in Fig. [2] 
and Fig. [3] are close enough to the relevant sulphur line that they 
would be indistinguishable. 

Our sample stars all have stellar parameters available in the 
literature (7^, logg, [Fe/H], and £ m i cro ). They are chosen to be 
cool giants in order to maximize the strength of the [S i] line. 
Since we are interested in halo stars with [Fe/H] < — 1 the un- 
resolvable Fe i blends expected to most seriously affect the 1045 
nm sulphur lines are small compared to the relevant sulphur lines 
and can safely be ignored. Our lines can thus be considered blend 
free and the measured equivalent widths are that of the sulphur 
lines only. 

Note that, as can be seen in Figs. [T] - [3] one should use these 
sulphur lines with caution for certain stellar parameters; e.g. for a 
solar-metallicity dwarf with a temperature of 4500 K the strength 
of the blending Fe i line would be roughly equal to that of the 
1045.5 nm Si triplet line. Also, the neighboring Cri and [Si] 
lines might become unresolvable and therefore blend when ob- 
serving at lower resolution. 

2.2. Reduction of the spectra 

The data were originally reduced using the CRIRES pipeline 
version 1.5.0^] The pipeline makes use of and automatically han- 
dles dark frames, flat fields from a halogen lamp and observa- 
tions of a ThAr lamp for wavelength calibration. To subtract the 
sky background the telescope was nodded between two positions 
along the slit during the observations. The two frames are sub- 
tracted from each other producing two stellar spectra in verti- 



2 http://kurucz.harvard.ed u/ 

3 http://www.graal.univ-montp2.fr/hosted/plez/ 
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Fig. 1. Equivalent widths of the 1082 nm [Si] line and potential blends for different stellar parameters. The leftmost panel shows 
the sum of the three CN lines listed in Tablepl The filled contours mark steps of 10 mA and the dotted contours mark the steps of 5 
mA in between. 



cally different places on the detector arrays. The pixels are sub- 
sequently added column wise in the slit direction and extracted 
for each spectrum exposure by the pipeline and then added to 
produce the final spectrum of the star. 

Unfortunately there is an optical ghosj^jin certain settings 
of the spectrometer, clearly visible as broad vertical bands in 
our fiat fields (see Fig. 4(a) and 



4(b) 



in our science frames (see Fig. 



but also as smaller dots 
where a smaller version 
of the ghost shown in Fig. |4(a)|can be seen). The optical ghost 



4(c) 



is present in some higher order settings for the echelle grating 
and is a result of retro-reflection from the detector to the grating 
which is redirected onto the detector in a different order. 

The ghosts in the flat fields lie in such a location on the de- 
tector that it will affect one of the two nod positions (compare 
the flat fields in Fig . |4( a) | and |4(b ) | wi th the stellar spectra in Fig. 



4(c) I. This results in a spectral artifact in the two spectral set- 
tings, most severely affecting our analysis of the [S i] line; by 
chance the spectral artifact lies right on top of the [S i] line in the 
pipeline reduced spectra making it impossible to measure (see 



Fig. 4(d) I. Therefore the data covering the [S i] line had to be re- 
reduced. We did this using the CRIRES pipeline version 1.11 .0^] 
using the flat field from the triplet setting instead, which was 
taken close in time. The spectral feature from the ghost is then 
situated in another part of the spectra not affecting our analysis 
(see Fig. 4(e) i. We found that this way of rescuing the data was 



at 



5 D escribed in the CRIRES User M anual 
http://www.eso.org/sci/facilities/paranal/instruments/crires/doc/ 

b Described in the CRIRES Pipeline User Manual, Issue 1.4 at 
ftp://ftp.eso.org/pub/dfs/pipelines/crire/crire-pipeline-manual-L4.pdf 



better than removing the ghost feature by dividing by a telluric 
star or by reducing the spectra without flat fields, thereby not in- 
troducing the ghost, and then dividing the spectra with a B star 
to remove the "pixel-to-pixel" variations in the detector. The lat- 
ter method was used successfully by Nissen et al. ( 2007| l when 



analysing science-verification spectra from CRIRES without flat 
fields. Our measured S/N values (per pixel of the spectra) of the 
reduced spectra reveal that our method indeed produces satisfac- 
tory signal-to-noise as shown in Table [2] 

The most relevant parts of the reduced spectra can be seen 
in Fig. [5] We normalized the spe ctra and fitte d the continua by 
using the IRAI^Jtask continuum ( |Tody|l993 i. 



2.3. Equivalent widths 

The equivalent widths of the 1045 nm Si triplet and the 1082 
nm [S i] line were measured using the splot task in IRAF by fit- 
ting a Gaussian profile to the line, or, in some cases, by straight 
numerical integration. For two stars (HD26297 and HD29574) 
the equivalent widths for the forbidden line were measured by 
de-blending due to the neighboring telluric line. The measured 
values are listed in Table [4] Since our spectra have high S/N the 
uncertainties in the measured equivalent widths are mainly due 
to uncertainties in the continuum setting. The quoted uncertain- 
ties have been estimated by fitting the continuum to achieve a 



7 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 2. Equivalent widths of the 1045.5 nm S i triplet line and its potential blends for different stellar parameters. The filled contours 
mark steps of 10 mA and the dotted contours mark the steps of 5 mA in between. 



maximum and minimum possible value for the equivalent width. 
The uncertainties in the measured equivalent widths will in the 
end produce uncertainties in the [S/Fe] less than the uncertain- 
ties introduced due to uncertainties in the stellar parameters (see 
Sect 1331. 



3. Analysis 

We have constructed stellar atmosphere models in spheri- 
cal geometry for the observed stars using the MARCS code 
(Gustafsson et al. 2008) and the stellar parameters quoted in 
Table [6] In all models we use an [a/Fe] value of +0.4. The 
MARCS code computes hydrostatic model photospheres under 
the assumption of LTE (Local Thermodynamic Equilibrium), 
chemical equilibrium, homogeneity and the conservation of the 
total flux. We used these models and the Uppsala Eqwi code (ver. 
7.06) together with our measured equivalent widths (see Sect. 



2.3 i to find two sulphur abundances for every star; one mean 



value based on the S i 1045 nm triplet and one based on the [S i] 
1082 nm line. The observed and synthetic spectra with sulphur 
abundances determined in the way described can be seen in Fig. 
[5] The synthetic spectra are calculated with the Bsyn code (ver. 
7.09), which is based on the same routines as the MARCS code. 
The continuous opacities are taken from the model atmospheres. 
Subsequently the spectra were broadened by convolving with a 
rad-tan function with a FWHM given in Table [6| in order to fit 
the data taking into account the macroturbulence of the stellar 
atmosphere and instrumental broadening. 



3.1. Atomic data 

The atomic data for the line list used in the model calculations 
are taken from the VALD database (Kupka et al. 1999), except 
the data for the S i triplet and [S i] fine which are taken from 
other sources, see Table [5] The [S i] line is a Ml intercombina- 
tion transition between the triplet ground state and the first ex- 
cited singlet state. The 1045 nm triplet is three relatively highly 
exited transitions in the triplet system. 

We note that the log(g/)-value of the [Si] line is debated 
Caffau et al. 2011) , and the derived 



(Asplund et al. 2009 



[S/Fe] values would typically be 0.08 dex lower using the older 
log(g/)-value used in Caffau et al. (2011 1 (they use log(g/) 
-8.617). 



3.2. Stellar parameters 

The observed stars all have Hipparcos parallaxes, but unfortu- 
nately the uncertainties for most of our observed stars are very 
large, see Table [T] This, in combination with the narrow wave- 
length coverage of CRIRES makes it impossible to determine 
the stellar parameters based on our observations alone, so we 
are forced to use stellar parameters from the literature. The stel- 
lar parameters were taken from Fulbright & Johnson ( 2003) , 
|Gratton et ah] (|2000), and |Burris et al.| ( |2000) (in turn mostly 



taken from |Pilachowski et aL ( 1996)) and in case a star is listed 
in two or more of those references the newest is chosen. 



Fulbright & Johnson (2003) obtained T e g from colours, 
[Fe/H] from Fe n lines, £ m i C ro from forcing all Fe I line strengths 
to give the same iron abundance, and log g from the standard for- 
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Fig. 3. Equivalent widths of the 1045.7 nm S i triplet line and its potential blends for different stellar parameters. The two leftmost 
panels show the sum of the three FeH lines and the two CrH lines listed in Table[3] The filled contours mark steps of 10 mA and the 
dotted contours mark the steps of 5 mA in between. 



mula: 



log — = log — + 0.4 (M v + BC - 4.72) + 41og -5? (1) 

80 M Q T G 



determining absolute magnitudes from literature values ( Hanson 
et al. 1998 Anthony-Twarog & Twarog 1994} in combination 
with fitting to a 12 Gyr isochrone and, where applicable, by us- 
ing Hipparcos parallaxes. Gratton et al. ( 2000[ ) determined T e g 
from colours, logg from the standard formula (Equation [T) with 
absolute magnitudes from |Anthony- Twarog & T warog (1994), 
£mkro from elimination of any trend for abundances as derived 
from Fei lines with different strengths, and [Fe/H] from Fei 
and Fen lines. Bur ris et aT] (|2"000 ) determined [Fe/H] from Fen 
lines, but took the rest of the parameters from Pilachowski et al. 
(1996) who obtained T e g from colours, ^micro from an iterative 
process fitting iron and calcium lines, and log g from three inde- 
pendent measurements: the standard formula (Equation [T} with 
absolute magnitudes from the Stromgren c\ -index, by assuming 
that the star's position in the M92 color-magnitude diagram, and 
through a derived average relation between surface gravity and 
effective temperature for metal-poor giants. The adopted stellar 
parameters are listed in Table [6] where we also give the uncer- 
tainties from the references. We note the slight inconsistency in 
adopting [Fe/H] -values as determined by plane-paralell model 
photospheres while we use spherical models. 

In our observed wavelength region there are four Fe i-lines 
which we used to check the [Fe/H] -values, given the other three 
parameters from the references are correct. Our [Fe/H]-values 



indeed fall within the uncertainties stated in the references. As 
can be seen in Fig. 5]the Fe and Cr lines in the vicinity of the [S 1] 
line, do not fit well with the given parameters for the two stars 
HD26297 and HD2957 4 Changing the T eS and [Fe /H] within 
the uncertainty stated in Fulbright & Johnson (2003 ) will how- 
ever make them fit. Despite this mismatch we have chosen to use 
the parameters determined by Fulbright & Johnson (2003 ) since 
their analysis is based on more and better known iron lines. 

Typical uncertainties for the derived sulphur abundances due 
to typical uncertainties in the stellar parameters used are listed 
in Table [7] As can be seen the 1045 nm triplet is most sensitive 
to the temperature, whereas the [S 1] line, coming from a tran- 
sition within the ground configuration, is not affected as much 
by uncertainties in the effective temperature, but more to uncer- 
tainties in the surface gravity through the dependence of the H 
bound-free opacity on electron pressure, and the iron abundance. 
All the measured sulphur lines are weak and therefore not much 
affected by the microturbulence. 



3.3. Non-LTE effects 



According to Takeda et al. ( 2005 ) the non-LTE corrections for 
the 1045 nm triplet are large and negative, i.e. the LTE abun- 
dances should be diminished in order to correct for the effects. 
They have calculated non-LTE corrections for a grid of model 
atmospheres with T eS = 4500,5000,5500,6000,6500,7000 
K, logg = 1.0,2.0,3.0,4.0,5.0, and [Fe/H] 
-4.0, -3 .0, -2.0, - 1 .0, -0.5 , 0.0, +0.5 . Their calculations 
use hydrogen collision cross sections based on the classical 
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(a) The flat field from 17 jan 2008 of de- 
tector two (of four) for the 1082 nm- 
setting. The optical ghost is clearly vis- 
ible as a broad band. 



(b) The flat field from 17 jan 2008 of de- 
tector two (of four) for the 1045 nm- 
setting. The optical ghost is clearly vis- 
ible as a broad band in another location 
on the detector. 



(c) The nodding of the spectrometer re- 
sults in two observations of the same 
star. The raw frames are here added to 
show both in one picture. The upper 
observation is affected by the flat field 
ghost but not the lower. The ghost can 
still be seen as the dots in the upper 
left-hand corner. The vertical line is in 
this case a telluric emission line. 




n n r.. I I ..1 n n t , , ■ 

1080 1081 1082 1083 1084 1085 1086 1080 1081 1082 1083 1084 1085 1086 

Wavelength (nm) Wavelength (nm) 

(d) The data for HD36702 reduced with (e) The data for HD36702 re-reduced using a 

pipeline showing the ghost-feature at 1082 flat field from 1047 nm showing the ghost- 

nm. feature at 1083 nm. 

Fig. 4. The optical ghost of CRIRES and its effects on reduced spectra. 



Table 6. Adopted stellar parameters for the observed stars. 



Star 


Teg 






[Fe/H] 




Refs. 


t b 

S macro 




(K) 




(cgs) 




(km s- 1 ) 




(kms- 1 ) 


HD 13979 


5075 ± 


100 


1.90+0.5 


-2.26 ± 0.09 


1.30 + 0.5 


1 


7.0 


HD21581 


4900 ± 


103 


2.24 + 0.2 


-1.64 + 0.09 


1.45° 


2 


4.8 


HD23798 


4375 ± 


82 


1.12 + 0.2 


-2.03 + 0.13 


2.20° 


2 


7.3 


HD26297 


4350 ± 


138 


1.46 + 0.2 


-1.51 + 0.11 


1.65° 


2 


5.4 


HD29574 


4200 ± 


119 


0.78 ± 0.2 


-1.70 + 0.12 


2.00° 


2 


6.0 


HD36702 


4366 ± 


100 


0.95 ± 0.2 


-2.06 + 0.11 


1.65 + 0.15 


3 


6.0 


HD44007 


4975 + 


103 


2.24 + 0.2 


-1.65 + 0.23 


2.20° 


2 


4.5 


HD83212 


4533 + 


100 


1.45 + 0.2 


-1.40 + 0.13 


1.80 + 0.15 


3 


6.0 


HD85773 


4450 ± 


100 


1.10 + 0.5 


-2.36 ± 0.09 


2.10 + 0.5 


1 


6.5 


HD103545 


4725 ± 


100 


1.70 + 0.5 


-2.14 + 0.09 


1.30 + 0.5 


1 


5.5 



References. (l) |Burris et allpOOO) ; (2) |Fulbright & Johnson|p003) ; (3) |Gratton eFaTlpOOOl i 

" The reference does not clearly state the uncertainty, so we assume a generic uncertainty of ±0.5 km s -1 . The size of this uncertainty does not affect 
the derived sulphur abundance much, see Table|7] since the sulphur lines used in this work are weak and do not depend much on microturbulence. 
* The macroturbulence is not taken from the references, but determined by us to fit the synthetic spectra to the observed spectra in Fig.pl 



formula by|Drawin ( 1968 1969 1 in the version of Steenbock & 



Holweger ( 1984JI i.e. S h — 1 is assumed. The work of Takeda 
et al.| ( |2005| ) is, to our knowledge, the only available non-LTE 



corrections for the 1045 nm triplet. |Korotin ( 2009| ) has, however, 
also performed non-LTE calculations for sulphur, but for the 
multiplet around 675 nm, the 869 nm doublet, and the 923 nm 



triplet. The overlapping parts of the two works agree well for 
lower-temperature stars, but differ for temperatures above 6000 
K. 

The non-LTE correction used ( [Takeda et aLj|2005) for our 
observed stars is shown in Table [9] Note that we have not at- 
tempted to interpolate in the grid, but chosen the corrections 
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Fig. 5. The observed spectra in black and synthetic spectra based on the equivalent width analysis in red. The normalized spectra 
are shifted by multiples of 0.2 along the flux-axis in the figure. The downward spikes sometimes occurring, e.g. close to the Fe line 
in HD83212, are due to bad pixels. 



from the models best applicable to the stars in this study. We note cause the populations of the energy levels are determined by col- 

the slight inconsistency of using spherical model photospheres lisions rather than radiative transitions. 

for our giants while applying non-LTE corrections calculated in 

plane parallel geometry. Further investigations on the effects of 

non-LTE corrections due to spherical geometry would be useful. 3.4. 3D modeling 

We have used the set of 3D hydrodynamic model atmospheres of 

Regarding the 1082 nm [S i] line, being an inter-combination red giants by Collet et aT| ( |2007| [2009 ) to perform spectral line 

resonance line, LTE is expected to be a good approximation, be- formation calculations for the [S i] line and the triplet. These cal- 
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Table 3. Potential blending lines turning up in our grid of mod- Table 4. 
els. lines. 



Measured equivalent widths for the observed sulphur 



Element 


Wavelength 


/t exc 


logfef) 




(nm) (air) 


(eV) 




Blends for the 1082.1 176 nm [S i] line 


Cri 


1082.1658 


3.013 


-1.678 


Fei 


1082.0882 


9.760 


-1.470 


FeH 


1082.0683 


0.619 


-0.331 


12 C 12 C 


1082.0866 


0.690 


-0.888 


12 c 12 c 


1082.1215 


0.610 


-1.316 


,2 C' 4 N 


1082.0779 


1.908 


-0.845 


12 C 14 N 


1082.0952 


1.552 


-1.379 


12 C 14 N 


1082.1382 


2.143 


-1.319 


Blends for the 1045.5449 nm S i 


line 


Pi 


1045.5870 


8.080 


0.000 


Fei 


1045.5404 


5.390 


-0.905 


12 C 12 C 


1045.5504 


1.370 


-1.233 


12 c 12 c 


1045.5614 


0.340 


-1.699 


CrH 


1045.5303 


0.600 


-0.698 


TiO 


1045.5562 


1.240 


-0.006 


TiO 


1045.5583 


1.170 


-0.344 


Blends for the 1045.6757 nm S i 


line 


Pi 


1045.6836 


8.080 


-1.000 


Cai 


1045.7098 


4.740 


-1.741 


Cri 


1045.6267 


4.420 


-1.855 


Cri 


1045.6467 


6.800 


-2.608 


Mm 


1045.6346 


5.820 


-1.210 


Fei 


1045.6455 


6.200 


-1.479 


Fei 


1045.6667 


4.190 


-4.669 


Fei 


1045.6939 


5.540 


-1.561 


12 C 12 C 


1045.6636 


0.250 


-1.854 


12 c 12 c 


1045.6731 


1.370 


-1.242 


12 c 13 c 


1045.7016 


0.270 


-1.345 


12 c 12 c 


1045.7133 


0.620 


-1.662 


12 C 14 N 


1045.6394 


1.160 


-2.931 


13 C 14 N 


1045.6915 


1.160 


-1.947 


13 C 14 N 


1045.7044 


1.530 


-0.864 


12 C 14 N 


1045.7253 


1.760 


-1.346 


CrH 


1045.6629 


0.078 


-0.728 


CrH 


1045.6983 


0.440 


-0.740 


FeH 


1045.6292 


0.100 


-2.539 


FeH 


1045.6379 


0.710 


-2.576 


FeH 


1045.7003 


0.470 


-0.869 


TiO 


1045.6596 


1.390 


0.352 


TiO 


1045.6653 


1.270 


0.042 


TiO 


1045.7061 


1.210 


-0.393 


Blends for the 1045.9406 nm S i 


line 


Vi 


1045.9169 


3.070 


-1.719 


12 C 12 C 


1045.8992 


0.620 


-1.672 


12 c 12 c 


1045.9890 


0.410 


-1.253 


12 C 14 N 


1045.9240 


1.330 


-2.576 


FeH 


1045.9672 


0.220 


-3.360 


TiO 


1045.8933 


1.510 


0.051 


TiO 


1045.9259 


1.210 


-0.074 



culations were compared to analogous ID calculations based on 
MARCS models to estimate the 3D-1D LTE corrections to the 
sulphur abundance for giants with different stellar parameters 
in order to check the effects of our neglect of convective inho- 
mogeneities on our abundance values. The results are shown in 
Table 1 

The parameters of the 3D models are not exactly the same 
as our sample stars, but these calculations provide a rough esti- 
mate of the 3D- ID differences and some qualitative conclusions 
can be drawn: the corrections for the [S i] line seem to grow 
larger for higher temperatures, larger surface gravity and lower 



Star 


1045.5 nm 

(mA) 


1045.7 nm 

(mA) 


1045.9 nm 
(mA) 


1082.1 nm 

(mA) 


HD 13979 


18 


± 1 


4: 


t 1 


12: 


1 1 


< 1 


HD21581 


48 


± 1 


16 d 


:0.5 


32 ± 


0.5 


6 ± 1 


HD23798 


28 ± 


0.5 






18 ± 


0.5 


10 ± 0.5 


HD26297 


39 


± 1 


12 


± 1 


26: 


t 1 


16 ± 1 


HD29574 


41 


± 1 


12± 1 


28 ± 1 


19 ± 0.5 


HD36702 


27 ± 


0.5 


6 + 


0.5 


19 ± 


0.5 


13 ± 1 


HD44007 


46 ± 


0.5 


15 d 


: 0.5 


32 ± 


0.5 


4± 1 


HD83212 


65 


± 1 


25 d 


: 0.5 


47 ± 


0.5 


17 ± 1 


HD85773 


16 i 


0.5 






12 + 


0.5 


8 ±0.5 


HD 103545 


26 


± 2 






14: 


t2 


< 5 



Table 5. Atomic data for the relevant sulphur transitions. 



Element 


Wavelength 
(nm) (air) 


(eV) 


log(g/) 


Transition 


Refs. 


Si 


1045.5449 


6.860 


0.250 


3 S°- 3 P 2 


1 


Si 


1045.6757 


6.860 


-0.447 


3 S?- 3 P 


1 


Si 


1045.9406 


6.860 


0.0285 


3 S°- 3 P, 


1 


[Si] 


1082.1176 


0.000 


-8.704 


3 P 2 - 1 D 2 


2 



References. (1) Zer ne et aE] ( |1997| >; (2) |Froese Fischer & Tachiev| 
pOTl) 

Table 7. A representative example (HD23798) of the effects on 
the sulphur abundances due to uncertainties in the stellar param- 
eters. 



Uncertainty 


Ale 


g e(S)triplet 


A log e(S) [S i] 


ST ef f - 


+ 100 K 




-0.11 


+0.04 


6\ogg 


= +0.2 




+0.08 


+0.09 


<S[Fe/H] = +0.15 


+0.03 


+0.05 




= +0.5 




-0.03 


±0.00 


Table 8. 3D - ID corrections for lo| 


;e(S). 


Teff 


l°gg 


[Fe/H] 


3D-lD [S i] 


3D-lD tl iplet 


(K) 


(cgs) 








4717 


2.2 


-1.0 


-0.04 


+0.17 


5131 


2.2 


-1.0 


-0.09 


+0.08 


4732 


2.2 


-2.0 


-0.08 


+0.19 


5035 


2.2 


-2.0 


-0.15 


+0.20 


4858 


2.2 


-3.0 


-0.27 


+0.20 


5130 


2.2 


-3.0 


-0.33 


+0.25 


4550 


1.6 


-3.0 


-0.09 


+0.17 



[Fe/H]. For our sample stars the 3D corrections for the sulphur 
abundances as derived from the [S i] line all seem to be around or 
below -0. 1 dex. The corrections for the triplet, on the other hand, 
are rather constant for the different 3D models and all seem to 
be around +0.2 dex. 

Our 3D analysis for these sulphur lines is the first for giants, 
but |Caffau & Ludwig" ( 2007|l explored the 3D correction for the 
[S i] line in the Sun and |Caffau eTal] ( p0"07] [20T0| explored the 
3D corrections for the 1045 nm triplet in the Sun, Procyon and 
four other dwarfs. All three works used C0 5 BOLD 3D hydro- 



dynamical atmospheres for calculating the 3D effects. Caffau & 
Ludwig (2007) found that the negative 3D-corrections for the 
[S i] line gets more significant (up to roughly -0.2 dex) for lower 
[Fe/H] and |Caffau et al. (2010 1 found positive 3D-corrections 
for the 1045 nm triplet more or less canceling the non-LTE cor- 
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rections applied (~ +0.1 dex). It is of course difficult to com- 
pare our models to the models used in Caffau & Ludwig ( 2007 ) 
and Caffau et al. (2007] [2010 ) because of the large difference in 
stellar parameters, but it can be noted, however, that all works 
predict a negative correction for the [S i] line and a positive for 
the triplet. 

We have also estimated corrections to the Fe abundances 
taken from the literature (see Sect. |3.2[ ), and found values from 
-0.1 dex to -0.2 dex for Fe I and from +0.05 dex to +0.1 dex for 
Fell. 



4. Results 

4.1. Sulphur 

The measured sulphur abundances and the applied non-LTE cor- 
rections from Take da et al.| ( |2005| l can be seen in Table [9] In 
this table we also list the difference between the LTE and non- 
LTE sulphur abundances derived from the 1045 nm triplet and 
those derived from the 1082 nm [S i] line. The plot of [S/Fe] vs. 
[Fe/H] can be seen in Fig. [6] and [7] We have used a solar sul- 
phur abundance of log e(S) & = 7.12 (Asplund et al. 2009} and 
all data in Fig.[6]are put on this scale. A different choice of solar 
sulphur abundance would simply shift all data systematically by 
e.g., 0.04 dex using the value of Caffau et al. ( 2011| l. 

In the Tableland Figs.[6]-[7]we have ignored 3D-corrections 
since we only have a coarse grid of models available for stars 
with parameters like our stars, but if they were to be applied they 
would work in the direction of lowering the sulphur abundances 
as determined from the [S i] line and adding to the abundances 
from the 1045 triplet. 



Table 11. Measured equivalent widths for the observed silicon 
lines. 



4.2. Other alpha elements 

In our narrow wavelength range there are few lines suitable for 
determining abundances for other a elements, but there are three 
Si i lines that can be used. The atomic data are taken from the 
NIST database and can be seen in Table [TUl 



Table 10. Atomic data for the silicon transitions used. 



Element 


Wavelength 




logte/) 


Transition 


Refs. 




(nm) (air) 


(eV) 








Si i 


1037.1264 


4.930 


-0.712 


jpo.Jg 


1,2 


Sii 


1084.3858 


5.862 


0.220 


ip_iD° 


1,3 


Si i 


1088.2809 


5.984 


-0.646 


3 D- 3 F° 


1 



Re ferences. (l)|Radziemski & Andrew|(|1965fr; 



(2)|Kelleher & Podobedova (20081; (3) Na har & Pradhan| |l993 ) 



Equivalent widths for the Si i-lines used are found in Table 

1 1 For stars with no values given for the 1037.1 nm and 1084.4 
nm lines, no observations were made and for the two stars with 
no values given for the 1088.3 nm line, the line was covered by a 
telluric line. The derived silicon abundances are shown in Table 

12 and the plot of [Si/Fe] vs. [Fe/H] in Fig. [7] We are not aware 
of any non-LTE calculations for these Sii-lines. 



5. Discussion 

Until a decade ago, only few studies of sulphur abundances in 
Galactic halo stars (or, more correctly, for stars more metal-poor 
than [Fe/H]< -1), had been made, mainly due to the lack of 



Star 


1037.1 nm 


1084.4 nm 


1088.3 nm 




(mA) 


(mA) 


(mA) 


HD 13979 




22 ± 1 


~4 


HD21581 


83 ± 1 






HD23798 




58 ± 1 


12 ± 1 


HD26297 




86 ± 1 


24 ± 1 


HD29574 




88+ 1 


21 ± 1 


HD36702 




63 ± 1 


14 ±2 


HD44007 


84 ± 1 


89 ± 1 


27 ±2 


HD83212 




123 ± 1 


44 ±2 


HD85773 




39 ± 1 


<5 


HD 103545 


57 ± 2 


63 ± 1 




Table 12. Derived silicon abundances. 




Star 


loge(Si) 


[Si/Fef 




HD 13979 


5.29 


0.04 




HD21581 


6.06 


0.19 




HD23798 


5.62 


0.14 




HD26297 


6.18 


0.18 




HD29574 


6.02 


0.21 




HD36702 


5.71 


0.26 




HD44007 


6.12 


0.26 




HD83212 


6.53 


0.42 




HD85773 


5.28 


0.13 




HD 103545 


5.70 


0.33 





" Using a solar silicon abundance of log e(S ;)o 
[20091 . 



7.51 JAsplund et"al] 



suitable diagnostic lines. The weak doublet around 869 nm, a 
transition from a relatively high level at 7.9 eV, was used but 
these lines are increasingly difficult to measure as the metal- 
licity decreases. However their formation is close to LTE con- 
ditions, and are therefore a recommended diagnostic when de- 
tactable ( |Takeda et al.|2005| >. In 2004 Nissen et al. and Ryde & 
Lambert used the triplet lines around 923 nm (Xexc — 6.5 eV) 
which are strong enough for determining sulphur abundances at 
high precision in halo stars, also for the metal -poor ones. The 
923 nm triplet is heavily affected by telluric lines, but this can 
nicely be handled with an observation of a telluric standard star 
( |Ryde & Lambert|2004"t |Nissen et al~|[2 004). A severe problem 
with these stronger lines has been that not all optical CCDs reach 
up to the wavelengths required. Also, the lines are presumably 
not formed in LTE (Takeda et al. 2005 ). Recently, the diagnostic 
power of new sulphur lines in the near IR, but beyond the reach 
of normal CCDs, have been explored by means of spectrome- 



ters like CRIRES at VLT and near-IR detector arrays (i.e., Ryde 
|20"06]|Nissen et al.|2007) . These lines are the 1082 nm [S i] line 
and 1045 nm triplet explored and used in this work. The 1045 



nm triplet has non-LTE corrections calculated by |Takeda et al. 



(2005 ) and the [S i] line, being a forbidden line, is formed under 
close to LTE conditions, which is confirmed by non-LTE calcu- 
lations by Korotin (private communication). This makes analy- 
ses based on other sulphur lines than the 869 nm doublet or the 
1082 nm [S i] line rely on more or less uncertairj^ but necessary, 
non-LTE calculations. 

In this paper we have analyzed ten halo giants observed in 
the near-IR at high spectral resolution in order to determine the 



It should be noted that non-LTE calculations are subject to a number 
of more or less important uncertainties and assumptions. The non-LTE 
corrections are therefore also uncertain. 
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Table 9. Derived sulphur abundances. 



Star 




Triplet 




[Si] 


Aloge(S) 


triplet- [S I] 




loge(S)" 


NLTE* 


[S/Fe] NLTE c 


loge(S) 


[S/Fef 


LTE 


NLTE 


HD 13979 


5.13 ±0.03 


-0.15 


0.12 


< 5.09 


< 0.23 


> -0.04 


> -0.19 


HD21581 


6.03 ±0.05 


-0.18 


0.37 


6.02 


0.54 


0.01 


-0.17 


HD23798 


5.65 ±0.04 


-0.18 


0.38 


5.48 


0.39 


0.17 


-0.01 


HD26297 


6.20 ±0.08 


-0.16 


0.43 


6.04 


0.43 


0.16 


0.00 


HD29574 


6.08 ±0.09 


-0.16 


0.50 


5.71 


0.29 


0.37 


0.21 


HD36702 


5.57 ±0.09 


-0.16 


0.35 


5.53 


0.47 


0.04 


-0.12 


HD44007 


5.89 ±0.02 


-0.15 


0.27 


5.86 


0.39 


0.03 


-0.12 


HD83212 


6.44 ±0.07 


-0.28 


0.44 


6.15 


0.43 


0.29 


0.01 


HD85773 


5.25 ±0.05 


-0.18 


0.31 


5.31 


0.55 


-0.06 


-0.24 


HD103545 


5.46 ±0.06 


-0.08 


0.40 


< 5.52 


<0.54 


> -0.06 


> -0.14 



" The uncertainty shows the line-t o-line scatter betwee n the three lines in the triplet. 

* non-LTE corrections taken from Takeda et al. j2005) . 

r Using a solar sulphur abundance of loge(S)o = 7.12 ( |Asplund et al.[ 2009 ). 
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Fig. 6. [S/Fe] versus [Fe/H] for the values deduced from the [S i] line shown with red circles and the non-LTE corrected values 
deduced from the 1045 nm triplet shown with blue squares. Typical uncertainties based on quadratic addition of the values in Table 
7] are shown i n upper right hand corner. The plus signs are measurements from |Chen et al.| ( |2002|l, upwards pointing tri angles from 
Caffau et al. ( 2005|l, a sterixes from Caffau et al. ( 2010| l, downwards pointing triangles from Takeda & Takada-Hidai ( 201 l| l, and 



exes from |Spite et al.lpOTT] ). They are all shifted to loge(S) G = 7.12 ( [Asplund et al.|2009| >. 



sulphur abundances from the near-IR triplet at 1045 and the for- 
bidden line at 1082 nm. We are able, for the first time, to analyse 
the 1082 nm [S i] line for halo stars and show that it is indeed a 

useful diagnostic in giants down to at least [Fe/H] 2.3. This 

result is in agreement with our model predictions (see Fig. [TJ 
showing that it should be usable even for lower [Fe/H] for cool 
stars with low surface gravity. Our model predictions also show 



that in dwarfs and subgiants this line is only detectable down to 

[Fe/H] 1 with the present means. 

Our results from the [S i] line favor a flat trend for [S/Fe] as 
a function of [Fe/H] for halo stars (see Fig. [6] and [7]). Fitting a 
line to the measurements^Jgives [S/Fe] = -0.021 ■ [Fe/H] +0.39 

9 We used the least square linear fit implemented in 
the routine fitexy presented in NASA's IDL library at 



11 



H. Jonsson et al.: Sulphur abundances in halo giants from the [S i] line at 1082 nm and the S i triplet around 1045 nm 



1.0 
0.8 
— 0.6 

<D 
L_ 

0.2 
0.0 



[S/Fe] from triplet 



HD85773 

HD13979 

, HD 103545 

HD36702 
' HD2379B 



HD29574 

HD44007 

.HD21513 

HD2B297 

HD832 1 2 



1 • 

9 : 



-2.5 



-2.0 -1.5 
[Fe/H] 



[S/Fe] from [S l] 



HD85773 

HD13979 

HD 103545 
HD36702 
HD2379B 



• ! 



HD29574 

i HD44007 

i ,HD215ie 

, „ HD26297 

HD832 1 2 



-2.5 -2.0 -1.5 

[Fe/H] 



[Si/Fe] 



HD13979 

i HD103545 

HD36702 
HD23798 



HD44007 
HD21518 

HD26297 

HD83212 



• * 



-2.5 



-2.0 -1.5 
[Fe/H] 



Fig. 7. [S/Fe] and [Si/Fe] versus [Fe/H]. The left panel shows LTE sulphur abundances deduced from the 1045 nm triplet in black 
squares and the non-LTE corrected values in red circles. The middle panel shows the sulphur abundances deduced from the 1082 
nm [S i] line. The right panel shows the silicon abundance. 



and the simple mean gives [S/Fe] = 0.43 with a standard de- 
viation of cr = 0.11. Since the [S i] line is not believed to be 
affected by non-LTE effects, and 3D effects are smaller, we con- 
sider these measurements more robust than the measurements 
using the 1045 nm triplet. By just comparing the two leftmost 
panels in Fig. [7] it seems that the non-LTE corrected sulphur 
abundances from the 1045 nm triplet are similar, but perhaps 
somewhat lower, than the abundances from the 1082 nm [Si] 
line. We have conducted a Kolmogorov-Smirnov tesj^] check- 
ing the probability for all our non-LTE corrected [S/Fe] 1045 
nm triplet measurements to come from the [S/Fe] distribution 
as described by our [S i] line measurements. The Kolmogorov- 
Smirnov statistic for these two distributions is 0.40 and p-value 
for this is 0.3 lp] meaning that formally the two distributions are 
similar to a low degree and therefore pointing toward a possible 
problem with the non-LTE corrections applied to the 1045 nm 
triplet. When instead comparing the sulphur abundance distribu- 
tion from the [S i]-line with the LTE sulphur abundance distri- 
bution from the triplet the Kolmogorov-Smirnov statistic is as 
expected larger; 0.50 and the p-value 0.11. Thus the non-LTE 
corrections of Takeda et al. (2005 ) at least adjust the LTE distri- 



bution of sulphur abundances from the triplet to be more like our 
resulting abundances from the [S i]-line. 

In our analysis we have so far ignored 3D-corrections due 
to the lack of models available for stars with parameters like our 
stars. It is of course very hard to inter- and extrapolate our coarse 
3D-grid to give us information applicable to the stars observed, 
but it would seem that the positive 3D corrections are large 
enough (almost +0.2 dex) to cancel the applied non-LTE cor- 
rections making the sulphur abundance measured from the 1045 
nm triplet higher than the abundances from the 1082 nm [S i] 
line. This would make the differences between the two data sets 



Press et al. 



(19921 



http://idlastro.gsfc.nasa.gov/ The routine takes the uncertainties 
in both [S/Fe] and [Fe/H] into account. 

10 We used the routine kstwo from 1 

11 The Kolmogorov-Smirnov statistic gives the maximum deviation 
between the cumulative distribution function of the triplet sulphur abun- 
dances from the values deduced from the [S l]-line. The p-value takes 
values between and 1 giving the significance level of the K-S statistic. 
A large p-value would, in our case, show that the cumulative distribu- 
tion function of our non-LTE corrected [S/Fe] 1045 nm triplet measure- 
ments is similar to the cumulative distribution function of our [S/Fe] 
[S i] line measurements. 



seen in Fig. [7] and in the results of the Kolmogorov-Smirnov test 
larger again and this might indicate that the non-LTE corrections 
of Takeda et al. (2005) are too small, or that our estimation of the 
3D corrections is too large. From Fig.|6]it indeed seems like the 
agreement between the two diagnostics would be even better if 
the [S/Fe] deduced from the [S i] line were lower and the values 
deduced from the triplet were higher, which at least coincides 
with the direction of the 3D corrections. It is however expected 
that non-LTE corrections fully incorporating convective inhomo- 
geneities would be different than the ID non-LTE corrections 
used. This needs further study with the coupling between 3D 
and non-LTE done self-consistently. We note that when remov- 
ing the applied 3D-corrections from Caffau eFaT] ( |2010) three of 
their four measurements, using the 1045 nm triplet, (asterixes in 
Fig.|6]l overlap very well with our 1045 nm triplet measurements. 
This is interesting since the measurements were made using the 
same instrument, the same diagnostic, and the same non-LTE 
corrections, but |Caffau et al.| (|2010)used dwarfs instead of giants 
like we do. This might be another indication that the non-LTE 
corrections work well in ID. 

Using the [S i] line as a diagnostic for determining sulphur 
abundance in halo red giants is preferable due to its small non- 
LTE and 3D effects, but using only this line is problematic, not 
only as exemplified by our data where an instrumental ghost un- 
fortunately hit our record of this line, but also because it might be 
blended by unidentified stellar lines, or telluric lines. Fortunately 
we cannot find any blends and, concerning the ghost, we found 
a way to rescue the data and determine the sulphur abundance 
from this line. Because of this risk more lines should be used si- 
multaneously, preferably both the 923 nm and 1045 nm non-LTE 
corrected triplets. It should be noted that the non-LTE formation 
of these lines are still not well understood (Take da et al.|2005] >. 
Observing both the 923 nm triplet and the 1082 nm [Si] fine 
would provide the same test of the non-LTE corrections of the 
923 nm triplet as presented here for the 1045 nm triplet. 

We find that the uncertainties in [S/Fe] in this study will 
mainly be due to uncertainties in the stellar parameters (in our 
case ~ 0.1 dex), the continuum drawing, (in our case < 0.1 
dex) and due to non-LTE corrections for the triplet and effects 
from convective inhomogeneities. Given these uncertainties no 
significant abundance trend in [S/Fe] vs. [Fe/H] but a scatter 
around [S/Fe] ~ 0.4 can be inferred from our data. To deter- 
mine whether this scatter is real we investigated the correlation 
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between [S/Fe] triplet and [S/Fe] [Sl] . They would be correlated 
if the scatter was real, but such a trend would also result from 
errors in the [Fe/H] used. We cannot find any correlation, but 
only scatter from what is expected from the uncertainties, ex- 
cept for HD 13979. We note that if we use the newer parameters 
and iron abundance of |Wylie-de Boer et al. ( 2010| l (T e ff = 5050 
K, logg = 2.25, [Fe/H] = -2.48, £ miao = 2.35 km s' 1 ) we get 
the sulphur abundances of HD 13979 to more resemble the trend 
of the others suggesting that our used [Fe/H] from Burris et al. 
(2000) might not be correct. 

Caffau et al] ( |2005| l find ~ 20 % of their analyzed stars to 



have 'high' values of [S/Fe] meaning that our randomly chosen 
sample of 10 stars would with 90 % probability find at least one 
from this population, if real, but we do not. Concerning the zig- 
zag trend of Takeda & Takada-Hidai (201 1 ) our sample of stars 
unfortunately do not go low enough in iron abundance to confirm 
or whether their rise is real or not. Our sample of stars do not 
have any stars in common with the previous studies to compare 
abundances from these with our [S i] abundances. 

The three Si i lines available in our spectra give a flat silicon 
abundance trend with roughly equal behavior as the flat sulphur 
abundance trend in our observed stars, which is what is expected 
from their common nucleosynthetic origin. The sulphur abun- 
dance of HD 13979 is low compared to the others, but also the 
silicon abundance is low, meaning that it is possible that this star 
is low in abundance of a elements, but, as mentioned earlier, the 
stellar parameters used for this star might not be correct. If we 
use the newer parameters of Wylie-de Boer et al. ( 2010[ ) also the 
silicon abundance of HD 13979 resemble the trend of the others. 

As can be seen in Fig. [5] the Fe and Cr lines in the synthetic 
spectra of HD26297 and HD29574 do not fit the data. This is 
likely due to the uncertainties in the adopted stellar parameters 
and the ^/-values of the lines. As noted earlier we still have cho- 
sen to use the parameters listed in Fulbrig ht & Johnson] ( |2003| l 
since our analysis of [Fe/H] would be based on Fe-lines with 
lower quality. 

An increasing number of investigations, but not all, are gath- 
ering evidence toward a view where sulphur follows the abun- 
dance trends vs. metallicity in the halo expected for the a ele- 
ments, i.e. at a constant value to low metallicities. This have im- 
plications for models of supernova and hypernova yields and can 
constrain such models. For example the hypernova yield models 
by |Nakamura eTal] pOOT) used by |Israelian & Rebolo| ( [200T] | 
to explain their 'high' values of [S/Fe] are not confirmed in 
newer yield calculations by Kobayas hi et al.| ([2006 ), but they 
predict a flat [S /Fe] vs. [Fe/H] trend for halo stars in agreement 
with e.g. our results. However, the predictions of high [S/Fe] 

energ 

include hypernovae with energies up to 30 x lO 31 erg, but it can 
be noted that Nakamura et al. (2001 ) also predict enhanced Si/O 
and Ca/O ratios for the most energetic hypernovae, which are 
not observed, so it seems that such very energetic hypernovae do 
not play any role for the chemical evolution above [Fe/H] 3. 

6. Conclusion 

In the literature, there are several proposed trends for the evo- 
lution of sulphur in the halo phase of the Milky Way. This has 
implications for our understanding of the nucleosynthesis of sul- 
phur and the sites where this occurs. The three currently dis- 
cussed trends are the flat trend of e.g., Spite et al. (2011 1, the 
scatter trend of|Caffau et al.||2010), and the recent zig-zag trend 



Nakamura et al. (2001 ) refer to hypernovae with expl osions 
;rgy of 100 x 10^ erg, whereas |Kobayashi et al.|(|2006[) only 



of |Takeda & Takada-Hidai| ( |201 1 



In trying to minimize the effects of non-LTE in the diag- 
nostic we have used the 1082 nm [Si] line for which LTE is 
a good approximation and found a flat trend for [S/Fe] similar 
to that in Spit e et al.| ( |2011| l and the trends expected for the a- 
elements. We also determined the sulphur abundance from the 
1045 nm triplet to investigate its usefulness in giants, especially 
concerning its sensitivity to non-LTE. Our empirical differences 
can be compared with existing non-LTE corrections. We find that 
with the non-LTE corrections from Ta keda et al.| ( |2005| l for the 
1045 nm triplet, the derived sulphur abundance in most giants 
are closer to the abundance determined from the forbidden line, 
but the data hints towards these non-LTE corrections being too 
large when ID models are used. When considering 3D we con- 
clude that they on the other hand might be too small. This, how- 
ever, is based on mostly qualitative conclusions from a coarse 
3D model grid, and it might be that the discrepancy between the 
two data sets is due to exaggerated 3D effects and not due to 
faulty non-LTE corrections. 

Not only since it is formed under LTE conditions, but also as 
the 3D corrections are smaller for the [S i] line as compared to 
the 1045 nm triplet, the 1082 nm [S i] line would be the preferred 
diagnostic of the two used in this paper. Even though the two 
diagnostics do not produce the exact same trend, our sulphur 
abundance from the near-IR triplet, corrected for the calculated 
non-LTE effects, also corroborates the idea of sulphur behaving 
like a normal a element and we do not find any 'high' values for 
[S/Fe] like in the scatter trend of |Caffau et al.| ( |2010| ). Neither 
do we find any 'high' values for [S/Fe] like in the zig-zag trend 
of Takeda & Takada-Hidai (201 1 ), but unfortunately our sample 
of giants does not go to [Fe/H] low enough to sample the region 
where their values of [S /Fe] rise. 

We have demonstrated that both the 1045 nm triplet and the 
1082 nm [S i] line should be useful for investigations of the halo 
phase of the Milky Way. The [S i] line has the advantage of not 
being subject to uncertain non-LTE corrections, but it also has 
the drawback of being a single line. An increasing number of 
investigations point toward sulphur behaving as a traditional a 
element in the trends in [a/Fe] vs [Fe/H] plots. However, there 
still exist considerable claims for higher values as the metallic- 
ity decreases. To convincingly determine the situation for the a- 
element sulphur, further systematic and homogeneous investiga- 
tions are needed using as many diagnostics and with as sophis- 
ticated analyses as possible (including non-LTE and 3D model- 
ing). 
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